On the joint inversion of SGG and SST data from the GOCE mission

The computation of spherical harmonic coefficients of the Earth’s gravity field from satellite-to-satellite tracking (SST) data and satellite gravity gradiometry (SGG) data is considered. As long as the functional model related to SST data contains nuisance parameters (e.g. unknown initial state vectors), assembling of the corresponding normal matrix must be supplied with the back-substitution operation, so that the nuisance parameters are excluded from consideration. The traditional back-substitution algorithm, however, may result in large round-off errors. Hence an alternative approach, back-substitution at the level of the design matrix, is implemented. Both a stand-alone inversion of either type of data and a joint inversion of both types are considered. The conclusion drawn is that the joint inversion results in a much better model of the Earth’s gravity field than a standalone inversion. Furthermore, two numerical techniques for solving the joint system of normal equations are compared: (i) the Cholesky method based on an explicit computation of the normal matrix, and (ii) the pre-conditioned conjugate gradient method (PCCG), for which an explicit computation of the entire normal matrix is not needed. The comparison shows that the PCCG method is much faster than the Cholesky method.